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Adiabatic dynamics of conduction electrons in antiferromagnetic (AFM) materials with slowly 
varying spin texture is developed. Quite different from the ferromagnetic (FM) case, adiabaticity in 
AFM texture does not imply perfect alignment of conduction electron spins with background profile, 
instead, it introduces an internal dynamics between degenerate bands. As a result, the orbital motion 
of conduction electrons becomes spin-dependent and is affected by two emergent gauge fields: one 
of them is the non-Abelian version of what has been discovered in FM systems; the other leads to 
an anomalous velocity that has no FM counterpart. Two examples with experimental predictions 
are provided. 
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I. INTRODUCTION 

Interplay between current and magnetization is an es- 
sential issue underpinning the field of spintronics 1 , which 
consists of two reciprocal problems: control of current 
through magnetization with a known configuration, and 
its converse, i.e., control of magnetization dynamics via 
applied current. In ferromagnetic (FM) materials with 
slowly varying spin texture m(r,t) over space and time, 
these issues can be solved by assuming that conduction 
electron spins always follow the background texture pro- 
file, known as the adiabatic approximation 2 ' 3 . The mi- 
croscopic basis underlying adiabaticity is the strong ex- 
change coupling H = —J<j-m(r,t) between conduc- 
tion electron spins and local magnetic moments, through 
which spin mistracking with the background causes large 
energy penalty and becomes highly unfavorable 4 . 

Under adiabatic approximation, the current - magne- 
tization interaction is recast into an emergent electrody- 
namics, in which its reciprocal influence boils down to a 
simple electromagnetic problem. Specifically, by diago- 
nalizing the local exchange Hamiltonian via local unitary 
transformation, fictitious electric and magnetic fields are 
generated into the orbital dynamics 4 7 
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where 0(r,t) and (j)(r,t) are spherical angles specifying 
the direction of m. As a consequence, the influence of 
background texture is represented by an effective Lorentz 
force F = sh(E+r x B) exerted on conduction electrons, 
where s = +1(— 1) denotes spin- up (-down) bands. The 
electric and magnetic components of the Lorentz force 
are responsible for the spin motive force 4 ' 6 and the topo- 
logical Hall effect, 8 respectively. In turn, back-reaction 
of the Lorentz force provides an interpretation to the 
current-induced spin torque exerted on magnetic tex- 
ture. 9 11 In a formal language, the adiabaticity induces 
an effective gauge interaction & n t = J^A^, where cur- 
rent 2\i acquires gauge charge according to s = ±1, and 



= ^4 M (m, dm) is the effective electromagnetic poten- 
tial representing space-time dependence of the texture. 
Variation over the current SCi nt / Sj^ = yields the effec- 
tive Lorentz force; and varying over the magnetization 
SCi n t /Sm = produces the spin-transfer torque. Thus 
the current-magnetization interaction is reciprocal. 

However, the above picture apparently fails in anti- 
ferromagnetic (AFM) materials where neighboring mag- 
netic moments are antiparallel. Conduction electrons are 
not able to adjust spins with local moments that change 
orientation on atomic scale. Nevertheless, the staggered 
order parameter n = (Ma — Mb)/2M s can be slowly 
varying over space-time, where Ma and Mb are the al- 
ternating local moments and M s denotes their magni- 
tudes. A natural question is whether a slowly varying 
staggered order still renders adiabatic dynamics of con- 
duction electrons in some other sense. This is desired 
knowledge for studying spin transport in AFM materials, 
especially the quest for current-magnetization interaction 
as that for FM materials. 

In spite of recent theoretical 12 and experimental 13 
progress, this problem has never been addressed. But, at 
the same time, AFM materials are believed to be promis- 
ing candidates for new thrusts of spintronics, 14 partly 
due to their tiny anisotropy, robustness against external 
magnetic perturbations, and absence of demagnetization, 
which brings prevailing advantages for experimental con- 
trol. In this paper, we develop the effective electron dy- 
namics in a smooth AFM texture by applying the non- 
Abelian Berry phase theory 15 17 on energy bands that are 
doubly degenerate. The physics of adiabaticity in AFM 
materials is found to be an internal dynamics between 
degenerate bands which can be attributed to a SU(2) 
Berry curvature. When translating into spin dynamics, 
the adiabaticity no more indicates spin alignment with 
the background, but a totally new evolution principle 
[Eq. (8a)]. Aside from spin dynamics, the orbital motion 
of conduction electrons is coupled to two different gauge 
fields: one leads to the non-Abelian generalization of the 
effective Lorentz force; the other results in an anomalous 
velocity that is truly new and unique to AFM systems. 
With comparisons to FM materials, this paper provides 
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a general framework on how a given AFM texture affects 
the dynamics of conduction electrons. The other side 
of the story, i.e., back-reaction of current on the AFM 
texture, will appear in a forthcoming publication. 

The paper is organized as follows. In Sec. II, the 
general formalism is presented where iso-spin is intro- 
duced. In Sec. Ill, our central results (8) are derived, 
followed by discussions on non-Abelian Berry phase and 
monopole, spin and orbital dynamics of conduction elec- 
trons, and comparisons of AFM electron dynamics with 
its FM counterparts. In Sec. IV, two examples are pro- 
vided, and the paper is summarized in Sec. V. Mathe- 
matical derivations are included in the Appendixes. 



II. FORMALISM 

Consider an AFM system on a bipartite lattice with 
local magnetic moments labeled by alternating Ma and 
Mb- The spin of a conduction electron couples to the 
local moments by the exchange interaction J(M /M s ) - cr, 
where cr denotes the spin operator of the conduction elec- 
tron, and M flips sign on neighboring A and B sublattice 
sites. In spite of antiparallel of neighboring moments, the 
staggered order parameter n = (M a~ M ~b) / ! 2M s usually 
varies slowly over space and time, and we can treat it as a 
continuous function n(r,t). Accordingly, the conduction 
electron is described by a nearest-neighbor tight-binding 
Hamiltonian locally defined around n(r,t): 
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where j(k) — —t ^ s e lk '° is the hopping term with S 
connecting nearest neighboring A—B sites (we set h = 1). 
In general J can be negative, but we assume a positive 
J throughout this paper. 

The local band structure can be easily solved as ±e(k) 
with e{k) = y/J 2 + |7(fc)| 2 , and in the adiabatic limit 
we neglect transitions between e and —e. Each of the 
two bands are doubly degenerate, and without loss of 
generality we will focus on the lower band — e with the 
two sub-bands labeled by A and £?, whose wave functions 
are \i/j a ) = e ik - r \u a ) and \i/j b ) = e ik - r \u b ). The Bloch 
waves \u a ) = \A(k))\t (r,t)) and \u b ) = \B(k))\± (r,t)) 
maintain local periodicity around (r,t), where 
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are local spin wave functions with = 0(r,t) and (j) = 
<fr(r, t) being spherical angles specifying the orientation of 
n(r, t). The periodic parts are spinors in the pseudo-spin 
space furnished by the A — B sublattices, 
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FIG. 1: A schematic view of Bloch waves in the lower band. 
Sub-band A means a local spin up electron has a larger prob- 
ability on the A sites and a smaller probability on the B sites; 
sub-band B means the opposite case. They are degenerate in 
energy and their wave functions have a finite overlap depend- 
ing on the ratio of J/e. 



which exhibit opposite spatial patterns schematically il- 
lustrated in Fig. 1. While (^al^b) = due to the orthog- 
onality of local spin eigenstates, {A(k)\B(k)) does not 
vanish, and we define this overlap as 



d{k) = (A{k)\B{k))- 
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which is a key parameter in our theory and £ < 1. It 
reaches maximum at the Brillouin zone (BZ) center and 
vanishes at the BZ boundary. From Eq. (4) we know £(fe) 
is a system parameter determined by the band structure, 
and it is constant since the energy conservation e = 
requires £ = 0. If J tends to infinity, the overlap £(fe) 
will vanish and the two subbands will be effectively de- 
coupled, by which the system will become a simple com- 
bination of two independent FM subsystems. 

We adopt the semiclassical approach to construct ef- 
fective electron dynamics, 15 where an individual elec- 
tron is described by a wave packet \W) = j dkw(k — 
k c )[c a \^a} +Cb\ipb)], where / dkk\w(k - k c )\ 2 = k c gives 
the center of mass momentum, and (W|r|W) = r c is the 
center of mass position. The coefficients c a and c b re- 
flect relative contributions from the two sub-bands and 
\c a \ 2 + |o>| 2 = 1. Since now the band is degenerate, non- 
Abelian formalism 3,15 17 must be invoked (see Appendix 
A), where dynamics between the A and B subbands in- 
troduces an internal degree of freedom represented by the 
isospin vector 



{ci,c 2 ,c 3 } 

{2Re(c a c 6 *),-2Im(c a c 6 *),|c a | 2 



lo>| 2 }. (5) 



The electron dynamics is characterized by the equations 
of motion of the three parameters & c , r c , and C, which 
can be obtained from variational principles with the effec- 
tive Lagrangian £ = (W\(i-§- t - H)|W). 3 ' 15 The detailed 
derivations are presented in Appendix A, here we only 
write down the results, 

(6a) 
(6b) 
(6c) 
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where r M = (t, r c ), but = (0, k c ) has no temporal 
component. In Eqs. (6) the • and x denote scalar and 
cross products in the isospin vector space. The Berry 
curvatures fi are obtained from the gauge potentials A. 
defined on the A and B sub-bands, for instance, 



= *<«i|0uK) 
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where r is a vector of Pauli matrices representing the 
isospin and i,j run between a, b. Other components of 
the Berry curvatures are explained in Appendix A. 



III. ELECTRON DYNAMICS 

Equipped with Eqs. (6), we are now ready to derive the 
dynamics of an individual electron. Note that the isospin 
vector C itself is not gauge invariant in the sense that dif- 
ferent choice of spin wave functions results in different C. 
Therefore, in deriving electron dynamics we need to re- 
late C to the real spin defined by s = (W\cr\W) (in unit of 
|) which is a physical variable and fully gauge invariant. 
While detailed derivations are lengthy and sophisticated, 
which is left for Appendix B, the final results are quite 
simple and elegant: 



s = (l-e)(s-n)h, 
k = • (Vn x s), 

r = -d k e - i(s x n) • h <9 fe ln£, 
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where h = dtn + (r • V)n, and we have omitted sub- 
script c of r c and k c for convenience of following dis- 
cussions. Equations. (8) are the fundamental equations 
of motion of a conduction electron in an AFM material 
with slowly varying texture, which are represented by 
joint evolutions of three variables (s,fc,r). An essential 
feature distinguishing the AFM electron dynamics from 
its FM counterpart lies in Eq. (8a), from which we know 
that the real spin s does not follow the background order 
parameter n in the adiabatic limit. 

Spin dynamics. The motion of s can be decomposed 
into a superposition of two motions: one strictly follows 
n (for stationary C) and the other represents mistrack- 
ing with n (for dynamical C), where the latter originates 
from dynamics between the A and B sub-bands and is 
unique to AFM materials (see examples in Sec. IV). It is 
worth emphasizing that the mistracking between s and 
n has nothing to do with any non-adiabatic process, but 
is entirely due to the non-Abelian nature of the prob- 
lem. The overall spin evolution can be attributed to 
the accumulation of a SU(2) non-Abelian Berry phase 
Vex.p[— i J A^-rdr^] along the electron trajectory 17 . As 
compared with its U(l) counterpart in FM materials 3 11 , 
which can be regarded as the magnetic flux of a Dirac 
monopole located at the center of the sphere spanned 




FIG. 2: (Color online) Left panel: the isospin vector C (blue 
arrow) in the local frame: C = c\0 + C2<fi + C3n, where and 
<t> are spherical unit vectors. Right panel: In our particular 
gauge, C (blue) is coplanar with n and s (red). The tip of 
C moves on a unit sphere, whereas tip of s is constrained on 
the ellipsoid whose semi-major axis is n and semi-minor axis 
having length £. 



by m, the SU(2) Berry phase here can be related to a 
't Hooft-Polyakov monopole in the parameter space. De- 
tailed discussions on the monopole can be found in Ap- 
pendix C, and the key point here is that a recent proposal 
of artificial 't Hooft-Polyakov monopole 18 can be realized 
in our AFM texture systems. 

Since |t(r,t)) and |^(r,t)) form local spin bases with 
quantization axis being n(r,t), the isospin C can be pic- 
tured as a vector in the local frame moving with n(r,t). 
In Appendix B, we have shown that in the particular 
gauge marked by x = (see Eq. (B18)), 

C = ciO + c 2 + c 3 n = + s 20) + s 3 n, (9) 

where {0, 0, n} form local bases associated with the local 
order parameter n(r,t) (Fig. 2, left). Equation (9) indi- 
cates two important properties: (1) C is coplanar with n 
and s, which is specific to the particular gauge \ = 0; (2) 
while the isospin vector is constrained on the unit sphere 



c + 
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1, the physical spin satisfies 



■ + «§ = !, 



~r ^2 °3 — x 5 UiiC piiy^J^cii opui dciijidiicd — 
which constrains the tip of s on an prolate spheroid with 
semi-major axis being n(r,t) and semi-minor axis on its 
equator having length £ (Fig. 2, right). The latter prop- 
erty is gauge independent and can be justified directly 
from Eq. (8a) without using Eq. (9) (see Appendix D), 
which, when written in gauge invariant from, becomes 
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So 
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For arbitrary gauges with x 7^ 0, it is easy to show that 
we always have S3 = C3 and s\ + s\ = i 2 (c\ + c|), but 
the angles between s\^ and will be different, i.e., C 
will not be coplanar with n and s. 

Now, the physical picture of adiabatic spin evolution is 
clear: as the background order parameter n(r,i) moves 
slowly in space- time, the prolate spheroid moves with 
it. The motion of physical spin s is a superposition of 
the relative motion on the spheroid and the motion of 



the spheroid itself. The overall motion of s described by 
Eq. (8a) is purely geometrical as dt can be eliminated 
on both sides, as a result, a given path of n uniquely 
determines a path of s on the spheroid which is indepen- 
dent of the Hamiltonian. Associated with this geometric 
motion, a SU(2) Berry phase is accumulated by the elec- 
tron wave function, which can be regarded as the (non- 
Abelian) gauge flux of a 't Hooft-Polyakov monopole at 
the center of the unit sphere spanned by n (Table I). 

A further remark: it seems to be a surprise that the 
magnitude of s varies on the spheroid since £ < 1, but 
how can the physical spin have a nonconstant magni- 
tude? We answer this question by studying the reduced 
density matrix for the spin degree of freedom. It is a 2 x 2 
matrix and can be written as p s = |(1 + a • cr), thus the 
expectation value of physical spin is s = Tr [p s cr] = a. 
Now since s 2 < 1, thus a 2 < 1, and what follows is 
Trp 2 < Trp s , which suggests that the electron is effec- 
tively in a mixed spin state. This can be attributed to 
the entanglement of spin and sublattice degrees of free- 
dom, specifically, because S3 = C3, we are able to infer 
the spin projection along n by measuring the probabil- 
ity difference on neighboring A — B sites (vice versa). 
The entanglement provides us with partial information 
of spin orientation from the knowledge of sublattice, this 
destroys full coherence of the spin states. 

Orbital dynamics. In correspondence with the novel 
spin dynamics, the orbital dynamics of an individual elec- 
tron also becomes non-trivial. By substituting Eq. (8a) 
into (8b) we get (see also Appendix B), 



k = (1 - £ 2 )(s • n)(E + r x B), 



E = -sin0(d t 0V0- V0d*<; 
B = ~ sin 6(V6 x V0), 



(11) 
(12) 

(13) 



the E and B fields here are the same as their FM coun- 
terparts, where they are responsible for the spin motive 
force 6 and the topological Hall effect, 8 respectively. Also, 
as in FM systems, it is easy to check that Eqs. (12) and 
(13) satisfy the Faraday's relation VxE+^ = 0. 

However, quite different from the FM case, the gauge 
charge s n in Eq. (11) is not just a constant, but involves 
the internal dynamics. In other words, the orbital motion 
is accompanied by a time-dependent gauge charge which 
should be determined by solving the coupled equations 
Eqs. (8) all together. Moreover, the factor £ 2 results from 
the non-commutative term 2*4^ x Af v in Eq. (7b), it also 
reflects the coupling between spin and orbital dynamics. 
The parameter £ E (0,1) plays a key role here: in the 
£ — >• 1 limit, 1 — £ 2 vanishes thus from Eqs. (8a) and (11) 
we get null results s = and k = 0. In the other limit 
where £ — >• 0, the solution of Eq. (8a) reduces to s = ±n 
if initial condition is s(0) = ±n(0), and Eq. (11) reduces 
to the FM Lorentz force equation, by which the system 
loses the manifest non-Abelian feature and behaves as 
two decoupled FM sub- systems. It deserves attention 
that in real AFM materials, both A and B sub-bands 
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TABLE I: Comparison of effective electron dynamics in FM 
and AFM textures. In the FM case, spin dynamics is trivial; 
along closed path Y the electron acquires an U(l) Berry phase 
which is the magnetic flux of a Dirac monopole; A Lorentz 
force is resulted in the orbital motion. In the AFM case, spin 
dynamics is non-trivial due to the mixture of degenerate sub- 
bands through a SU(2) non-Abelian Berry phase, which is the 
gauge flux generated by a 't Hooft-Polyakov monopole; The 
orbital dynamics is subject to a spin-dependent Lorentz force 
and an anomalous velocity. 



host majority carriers, but they are subject to effective 
Lorentz forces of opposite directions, which may lead to 
non-trivial spin transport. 

Furthermore, the real space dynamics governed by 
Eq. (8c) exhibits spin-orbit coupling through the anoma- 
lous velocity term - (s x n) n dk In £. It is along the same 
direction as <9&£, so Eq. (8c) amounts to give a modified 
group velocity. We mention that this term is unique to 
AFM textures and has nothing to do with the anomalous 
velocity studied in FM or quantum Hall systems 3 . It orig- 
inates from the Qjjf v curvature that joints real space with 
BZ, the importance of which has been overlooked before. 
For better comparison, we summarize the fundamental 
electron dynamics of FM and AFM textures in Table I. 

A final remark on the theory part: In real materials 
with impurities, our fundamental equations (8) are valid 
so long as spin coherence length is as large as, if not 
more than, the typical width of the texture. While this 
is quite true in FM materials, its validity in AFM mate- 
rials awaits experimental verification. At extremely low 
temperatures, spin-flip scattering is dominated by mag- 
netic impurities which can be made negligibly small in 
clean samples. Besides, spin-independent scattering pro- 
cesses (e.g., electron-phonon scattering) do not destroy 
our essential conclusions if r is understood as the drift 
velocity of carriers. We mention that AFM spintronics 
is an emerging field where very little is known. While 
it shares some similarities with the established FM spin- 
tronics, it is not always correct to copy ideas from FM 
systems. The adiabatic electron dynamics studied in this 
paper is one example. 



IV. IV EXAMPLES 

First, consider a spiraling AFM texture sandwiched 
by two ferromagnetic layers (see Fig. 3). This magnetic 
structure has been realized in Co/FeMn/Py trilayers in 
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FIG. 3: (Color online) Left: FM/AFM/FM trilayer with op- 
posite FM orientations on two sides. The black double ar- 
rows represent the A-B sublattices of the AFM layer, which 
is dragged into a spiraling texture due to exchange bias on the 
interfaces. Right: incoming electrons only enter the A sub- 
band due to the upper FM polarizer, the out-going electrons 
partially occupy the B sub-band depending on the value of £. 

a recent experiment 19 , where the FM order of Co layer 
is nearly fixed but that of Py can be rotated by external 
magnetic field. The AFM order is dragged into a spiral 
due to the exchange bias effect on the AFM/FM inter- 
faces. The layer thickness of FeMn is roughly 10 ~ 20 nm 
and can be made even larger, which far exceeds the lattice 
constant thus adiabatic approximation is valid; Mean- 
while, typical spin coherence length is larger than the 
layer thickness at low temperatures so that spin evolu- 
tion is governed by Eq. (8a). 

When an electron flows from top to bottom with ap- 
plied current, the top FM layer polarizes its spin so that 
it enters the A sub-band across the interface. According 
to Eq. (8a), the physical spin orientation of the electron 
after passing through the AFM layer is rotated by II = 
7r — arctan[£ tan£7r] if £ < ^, and II = — arctan[£ tan£7r] 
is £ > ~. This is a topological result that only depends 
on the initial and final directions of n, but is independent 
of the texture's profile detail. When £ — ^ 0, II reduces to 
7r, which means the electron spin follows n and remains 
in the A sub-band, thus it flows into the bottom FM layer 
with a lower resistance; in the £ — >• 1 limit, II vanishes 
and the electron completely evolves into the B sub-band 
thus experiencing a higher resistance. For an arbitrary £ 
and an arbitrary total rotation of the spiral denoted by 
<£, the electron will partially evolve into the B sub-band 
with the wave function cos(£<I>/2)|?/; a ) + i sin(£<I>/2)|^}, 
thus the total resistance is 

P = Po + -Ap[l-cos(£$)], (14) 

where po is the intrinsic resistance of the AFM texture 
itself, which depends monotonously but not too much on 
<£. Ap represents the magnetoresistance of the spin valve 
which is determined by material details of the two FM 
layers and is independent of <£>. If $ is increased beyond 
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FIG. 4: (Color online) Spin evolutions for three different £'s 
when n(t) is moving round a cone with constant angle 6 from 
the z axis. Upper panels: the tip of s respects two constraints: 
it stays both on the cone's bottom (small gray slab) and on 
the spheroid described by Eq. (10) (blue ellipsoid), thus the 
vector s is confined in between two cones with different semi- 
angles. Lower panels: orbits of the tip from bird's eye view. 
The topology of the orbits is separated into two classes (left 
and right) by the critical case (middle) where the inner cone's 
semiangle shrinks to zero. The orbits are not necessarily com- 
mensurate with n. 

7r, p will reach a maximum at <& m = 7r/£ and then re- 
duces. The resistance maximum, if observed, serves as 
an experimental verification of Eq. (8a). Moreover, mea- 
suring <I> m also enables us to find £ without calculating 
the band structure. 

We remark that the above results survive in the pres- 
ence of diffusive processes so long as spin-flip scattering 
is ignored. The reason is that spin-independent scatter- 
ing only deflects fc-space orbit, whereas the s dynamics 
is determined by the variation of n that is blind to k in 
one dimension. In addition, FeMn is a non-collinear anti- 
ferromagnet that has more than two sub-lattices. To test 
our theory unambiguously, we can replace FeMn by the 
collinear IrMn which is feasible for current technique 20 . 
Moreover, we are aware of the experiment 21 where the 
spiraling AFM texture exhibits spatial periodic patterns, 
it provides a better way of realizing large &'s. 

Consider a second example where n(t) is varying round 
a cone of constant semiangle in the laboratory frame, 
which can be realized in a spin wave (see Fig. 4). Accord- 
ing to Eq. (8a), we know that ds z = due to dn z = 0, 
thus, the tip of s should stay in the bottom plane of the 
cone. On the other hand, we learn from Eq. (10) that 
the tip is constrained on the spheroid that moves with 
the instantaneous n{t). Therefore, the actual orbit tra- 
versed by the tip is contained in the intersection of the 
two constraints. Through some straightforward geomet- 
ric analysis, we know that s is bounded between the n 
cone and an inner cone whose semiangle depends on £. 
Figure 4 depicts the actual orbits of s for three different 
£'s: they all exhibit precession and nutation, which can 
be easily read out from the bird's eye view. Remark- 
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ably, the motion of s falls into two topologically distinct 
classes separated by the critical condition 



Sc 



cos 2 



(15) 



In a real spin wave, is nearly zero, thus £ c « l/y/2. For 
real materials, we expect t < J, thus from Eq. (4) we 
know that for a partially filled band, £ is always smaller 
than l/y/2. Therefore, the £ < £ c phase is more realistic. 



V. CONCLUSIONS 

In this paper, we find that a slowly varying AFM tex- 
ture renders adiabatic dynamics of conduction electrons, 
which are described by three coupled equations of motion 
[Eqs. (8)]. Quite different from the FM case, the adia- 
baticity in AFM materials does not imply strict align- 
ment between conduction electron spins and the pro- 
file of background texture. Instead, the adiabatic spin 
evolution is a superposition of a motion following the 
background order plus a motion on a prolate spheroid 
attached to the local order, where the latter originates 
from internal dynamics between degenerate bands. The 
overall motion of the spin is still geometric; it can be 
attributed to the accumulation of a SU(2) non-Abelian 
Berry phase originating from the gauge flux of an effec- 
tive 't Hooft-Polyakov monopole in the parameter space. 

The corresponding orbital dynamics shares some simi- 
larities with FM materials in that the fc-space dynamics 
can be described by an effective Lorentz force equation. 
However, two prominent differences in the orbital dy- 
namics distinguish an AFM system from its FM coun- 
terpart: first, the gauge charge is dynamical rather than 
constant, by which spin and orbital motions no longer 
separate; second, the group velocity is renormalized by a 
spin-dependent anomalous velocity, which is quite differ- 
ent from what has been studied before. 

Our theory lays the foundation for charge and spin 
transports in AFM texture systems, which will be ap- 
plied to real materials in the future. The validity of our 
theory needs to be tested experimentally since available 
data on AFM spintronic materials are very rare. Theo- 
retically, this paper solves only the first half of the whole 
story; the other half, i.e., the converse effect regarding 
the back-reaction of current on background AFM order, 
will appear in a forthcoming paper. 
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Appendix A 



For simplicity, we denote the joint space-time coordi- 
nate as r M = (£, r c ), so that the local Hamiltonian H(r M ) 
is parameterized by r M and the wave packet of the elec- 
tron can be expressed as 



\W{r^k c ))=J dkw{k - k c )[c a \ijj a ) + c 5 |^ 6 )] 

= J dkw(k - fe c )[ca|A(fc))|t(r M )) + c b \B(k))\Ur^))}, 

(Al) 

where the profile function w(k — k c ) satisfies J dk\w(k — 
k c )\ 2 = 1 and j dkk\w(k — k c )\ 2 = k c . The coefficients 
c a and C6 reflect relative contributions from the two sub- 
bands. Here, we define the isospinor c = [c a ,c b ] T , then 
the normalization condition \c a \ 2 + |q,| 2 = 1 becomes 
c*c = l. 22 The effective Lagrangian for the wave packet 
is then written as 15 

C=(W\{i^-H)\W) = £(r li ,k lt ,c-, f„,k^5) 



where fc M 
trast to r u 



e + k c - r c + iPc + c^A^r^ + A£fc M )c, (A2) 

= (0, fe c ) has no temporal component in con- 
(t, r c ), but it is still written this way just to 



simplify symbols. Here the Berry connections are 2x2 
matrices and are functions of both r M and . The real 
space components are defined by 15 17 



(u a \d^\u a ), {u a \&l\u h ) 
{ u b\dl\u a ), {u h \dl\u h ) 



cos 0d L1 



i{-id^0 - sin 0d^ 



-Ti£sin 0d^ 



sin 0d^ (j)) — cos 0d^ <; 



■T2f0 M + T 3 COS0d M 0], (A3) 



where {ti,T2,T3} = r is a vector of Pauli matrices rep- 
resenting isospin, it should not be confused with the real 
spin operator a. Similarly, the momentum space compo- 
nents are 



(u>a\d*\u a ), (u a \d*\u b ) 
(v>b\d*\u a ), (u b \d*\u b ) 














(A4) 



which in general should not vanish if spin-orbit coupling 
terms are added to the original Hamiltonian. To avoid 
complicated manipulations on matrix products in the fi- 
nal results, we define the isospin vector 



C = {ci,c 2 ,c 3 } = cVc 

= {2Re(c a c 6 *),-2Im(c a c 6 *),|c a | 2 



|Q>| 2 }, (A5) 
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with which the Berry connection A 1 ^ can be expressed as 
a vector in the isospin vector space (the adjoint represen- 
tation), 



Al = -Tr[rAl] = {(Al)\ (^) 2 , (A^f) 



{-£ sin0d M 0, £8^0, cosOd^}, (A6) 



where we have used |Tr [r^Tj] = 5ij. In a similar way, 
A 1 ^ = ^Tr[rA^] = 0, but we still keep it at this step for 
the completeness of the entire formalism. The effective 
Lagrangian now becomes 

C = e + + ic^c + C ■ [A^ + A^k^ (A7) 

where the ic^c term cannot be expressed in terms of C 
and C, but this poses no problem in the following. Both 
Eqs. (A2) and (A7) are useful. 

To obtain the equations of motion for the three pa- 
rameters (C, r M , k/j), we resort to the following variational 
principles. 

(1) SC/Sc = on Eq. (A2) gives: 



(A8a) 
(A8b) 



c — A^yTyL + A^k^c 
gt = _ict[A^ + 4^] 

then from C = c^tc we have, 

C = cVc + C^TC 

= ir^[rAl - A^r)c + ^[r^ - Ajr]c (A9) 

where ri M terms are matrix products. In view of the 
decomposition Eq. (A6), we take a specified component 
of Eq. (A9), 

Ca = ir^A^C^TaTp - Tf3T a )c 

+ ik^A^Yc^TuTp - T/3T a )c 

= -2e afh [{A r J% + (A*)%](<*Trf, (A10) 

when written in the isospin vector form, the above equa- 
tion becomes 



C = 2Cx{A r ^ + A%) 

which proves Eq. (6a). 

(2) SC/Sr^ = on Eq. (A7) requires some care: 



(All) 



5C_ _ dC_ _ d DC 

<5r M dr^ At dr^ 



with, 



= d;s + C ■ [(d*AZ)r v + (d;At)K\ (A12a) 

= fc M + 2[r„C • (Al x ^) + KC ■ (At x ^)] 
+ C-[{dlAl)r v + {d k v Al)L] (A12b) 



where in the last line Eq. (All) has been used. Hence 
we obtain the equation of motion 

^ = ^e + C-[n^r„ + n^], (A13) 

where the Berry curvatures are defined as 

= d^A: - dlAl + 2AI x Al, (Al4a) 

n r Jl^d;At-dtA^ + 2A^xAt, (A14b) 

they are antisymmetric tensors with permutations of r M 
and and at the same time they are vectors in the 
isospin vector space - the internal space unique to non- 
Abelian gauge theory. 

(3) bC/dk^ = on Eq. (A7) follows quite similar pro- 
cedures as above, and the equation of motion is: 

r li = -d*e-C-[n%,r v + n%k v ], (A15) 
where the Berry curvatures are 

n% = - d k v Al + 2AI x Al, (Al6a) 

= d'Al - dlA\ + 2A\ x Al, (A16b) 

which, together with Eqs. (A14), form a generalized ma- 
trix of Berry curvature jointing real space and BZ into a 
unified parameter space, 



C^kr C^kk 



(A17) 



Eqs. (A13), (A15), and Eq. (All) justify the main conclu- 
sions in Sec. II. It is worth mentioning that the Abelian 
version of the above formalism is a well-established field, 
whose great success has been proved by compelling ev- 
idences in the last decade. 3 However, the non- Abelian 
version is still underdeveloped in previous work. 15 



Appendix B 

Before deriving Eqs. (8), special attention should be 
paid on the fact that gauge fields (Berry curvatures) in 
non- Abelian theory are not gauge invariant, but gauge 
covariant. It is the isospin scalar s C • appearing in 
Eqs. (6) that respect gauge invariance. Specifically, as 
we make a gauge transformation on the wave functions 
\ip a ) and l?/^), change of ft just compensates that of C. 
However, to perform real calculations, we have to work 
in a specified gauge. In this paper, the gauge is fixed in 
the choice of local spin eigenstates, which are obtained by 
acting U(r,t) = e -^^/^ e -^y0/2 e -ia z x/2 on the eigen _ 

states of a z . While 6(r, t) and t) are physical, x( r : i) 
is not and can be chosen arbitrarily; the gauge is fixed 
by setting x — 0. 

From Eq. (A6) the real-space curvature is obtained: 

= ®\iA v — d r v A^ + 2*4^ x A v 

= {0, 0, (e - 1) 1 - smeid^d^ - died^)} 

= {0, 0, (e 2 -l)in-(^nx5>)}, (Bl) 
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where n = {sin cos 0, sin sin 0, cos 0} is the local or- 
der parameter. We see that only the third component is 
non-zero in our particular gauge marked by x = 0- But 
one can check that in any gauge with x / 0, the first 
two components do not vanish. However, the third com- 
ponent is actually gauge invariant and it has the form of 
skyrmion density. The cross components of Berry curva- 
ture are obtained in a similar way, 



dtA:,+2A: t xAi 



-{dUsmed^ -d^o;e, o}, (B2) 



where again the first two components are changeable sub- 
ject to gauge transformations, whereas the third is gauge 
invariant. Moreover, due to = 0, the BZ space Berry 
curvature vanishes and we will not mention it in the 
following. Now let us turn to the equations of motion. 

(1) Substitute Eq. (A6) into Eq. (All), noting that 
6 = r^d^O = d t + r c -V0 and the same for 0, 



_d 

dt 



Our target is to transform the dynamics of C to the dy- 
namics of the physical spin defined by 







cos 00 -£0 




Cl 




C2 




— cos 00 — £ sin 00 




C2 


. (B3) 


C3_ 




£0 £ sin 00 




_C3_ 





(B4) 



which respects gauge invariance. From Eqs. (2), (Al), 
and (A5), and (B4) we know the components of s in the 
lab frame after some tedious algebra, 



x = cs sin cos + £[ci cos cos — C2 sin 0] , (B5a) 
y = cs sin sin + £[ci cos sin + C2 cos 0] , (B5b) 



C3 cos0 — £ci sin0, 



(B5c) 



then we take the total time derivative over each compo- 
nent of s, for example, 



Sx — cs sin cos + cs (cos cos 00 — sin sin 

+ £[ci cos cos — ci (sin cos 00 + cos sin 

— £[C2 Sin + C2 cos 00] 

= 03(1 — £ 2 ) (cos cos 00 — sin sin 
= c 3 (l -£ 2 )r^ 



(B6) 



where in deriving the second equality above Eq. (B3) has 
been used. Similarly, 



s y = c 3 (l-£ 2 )r^, 



c 3 (l-e)n z . 



(B7) 



To eliminate C3 in the above equations, we reverse 
Eqs. (B5) and obtain, 

C3 = s x sin cos + s y sin sin + s z cos 6 = s ■ n (B8) 



then from Eqs. (B6) and (B7), we obtain a simple and 
elegant equation of motion for the physical spin, 



s = (l-e)(s-n)h, 



(B9) 



which justifies Eq. (8a). 

(2) To justify Eq. (8b), we substitute Eqs. (Bl) and 
(B2) into Eq. (A13), regarding that fc M has only spatial 
but no temporal components, we arrive at, 

k c = 7^3 (£ 2 - l)n • (Vn x h) + -£[c x sin0V0 - c 2 V0] 
1 



= -n-{Vnx[(^-l)( S -n)n]} 
= — • (Vn x s) 



(BIO) 



where £ = and Eq. (B9) have been used. We also have 
ignored d^e term in Eq. (BIO) since the band structure 
is only a function of k and independent of space-time in 
the adiabatic approximation. As in Sec. Ill, we will omit 
the superscript c in the following for simplicity. To make 
better comparisons with the spin motive force 4 ' 6 and the 
topological Hall effect 8 discovered in FM materials, we 
also derive another suggestive form of Eq. (8b). Take an 
arbitrary component i of Eq. (BIO), 



^ = lc 3 (£ 2 -l) sin 0(5,1 



-c 3 (e-l)sm0{[d z 0d^-d t 0d^} 

\cz(t? - l){sin0[d,0d t - d t 6d^\ 
+ sin OsijkSkimrjdiOdmcj)} 



(Bll) 



where = dtO + r^0 (the same for 0) and the identity 
SijkSkim = SuSjm-dimSji have been used. Eq. (Bll) can 
be written in a concise way as, 

k=(l-£ 2 )(s-n)(E + rxB), 



E = -sin0(<9 t 0V0- V0<9^ 



B 



sin0(V0 x V0), 



(B12) 



(B13) 



which proves Eqs. (11) - (13) in Sec. III. 
(3) From Eq. (A15) and Eq. (B2), we have, 

r = -d k e + ^<9fc£(ci sin 00 - c 2 0). (B14) 

There is a smart way to eliminate in terms of physical 
spin. Notice that in the special gauge marked by x — 0, 
the isospin vector C can be pictured as a vector in the 
local frame extended by 0, 0, and n (See Fig. 2): C = 
c\0 + C20 + C3n, and components of the real spin s in 
this local frame are 

s\ = s x cos cos + Sy cos sin — sin = £ci (B15) 

52 = — s x sin0 + s y cos0 = £c2 (B16) 

53 = s x sin cos + sin sin + s 2 cos = C3 (B17) 
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where Eqs. (B5) have been used, and we obtain the im- 
portant relation, 



s = £(ci0 + c 2 0) + c 3 n. 
Since h = 00 + sin we immediately have 

n ■ (s x ft) = £(ci sin 00 — C20), 
thus Eq. (B14) becomes 

r -dkS - ^(s x n) • n<9fcln£ 
which justifies Eq. (8c). Moreover, 

1 - e d k e 



<9 fc ln£ 



(B18) 



(B19) 



(B20) 



(B21) 



thus the term ~(s x n) • n<9fc ln£ has the same direction 
as the group velocity, it represents the modification of 
group velocity due to AFM texture. 

A further point should be added is that in the most 
general case, the effective Lagrangian Eq. (A2) should 
also contain a term representing self-rotation of the wave 
packet —lm[ci(dJ l Ui\(e — T~0\d'^Uj)cj] 1 but after some so- 
phisticated manipulations one can show that this term 
vanishes for similar reasons as the vanishing f2j^. 

Appendix C 

To study the monopole, we should turn to a differ- 
ent coordinate system. By assigning a variable magni- 
tude to n, we define the dimensionless order parameter 
R = Rn = jfi. Then the Berry connection can be 
equivalent ly defined in the R space, which relates to the 
original one by A^dr^ = AidRi, where 



An i 



CM ait> e<tiau> 
eaiAit) aiAii) 



(Cl) 



where £ = |7|/ ^/ R? + I7I 2 is also a function of i?, and 



l7l=Ei 



ikd 



depends on the position in BZ. Written in 



spherical coordinates, components of Eq. (Cl) are 



A 



n 



O,A 



2R 



-i 

1 



,Aa 



1 

2# 



cot# 



-e 

- cot ^ 



(C2) 

To see the monopole, we should further make a singular 
gauge transformation on the potential 23 , 



A' e = SA e S^ + i\-SdgS^ 
R 

2i? 

a; = s^s* + » 

2i? 



1 



i?sin6> 



e -i<p 
b 



— sin 6 e l ^ cos 
e 1 ^ cos # sin 9 



(C3a) 



(C3b) 



with the unitary matrix 

,/2. 

5 = 



-e"^ /2 sin • 
e^/ 2 cost 



(C4) 



Finally, expressing the gauge potential in Cartesian co- 
ordinates, we obtain: 



A x = A'n cos cos 6 — A' * sin 6 = 



(i-O 

2i? 2 



y iz 

-iz -y 



A y = A' e sin (j) + cos < 



A. = An sin 



(1-0 
2i? 2 



(i-O 

2R 2 


-2/ + ix 



-x z 

Z X 

-y — ix 




they are nothing but the "hedgehog" gauge potential of 
a 't Hooft -Polyakov monopole at R = 23 , 



A' = 



i - OT 

2R 2 



^ijkRj&ki 



(C5) 



where Ri = {x,y,z} and cr^'s are Pauli matrices. The 
radial profile of Eq. (C5) is determined by the factor 1 — 
£(R). For fixed nonzero 7, 1 — £(i2) tends to 1 as R —> 00, 
and 1 — ~ i? 2 as — >• 0, which cancels the i? 2 in the 
denominator thus the gauge potential is regular at origin. 
The form and behaviors of Eq. (C5) are all the same as 
the gauge potential proposed by Sonner and Tong for 
realizing artificial 't Hooft-Polyakov monopole 18 . In fact, 
our parameter £ can be understood as the f(B) factor in 
Ref. [18], they both reflect the overlap of (partial) wave 
functions from doubly degenerate bands. However, at BZ 
boundary 7 is zero, thus £ = regardless of R. In this 
case, Eq. (C5) becomes the gauge potential of a Wu-Yang 
monopole and the origin R = becomes singular. 

Along with the gauge potential, we can define the as- 
sociated Higgs field, 



'(A\a 3 \A) (A\a 3 \B) 
(B\a 3 \A) (B\a 3 \B) 



= vr^e * 3 (ce) 



which is the pseudo-spin polarization of AFM systems. 
It describes the extent to which the conduction electrons 
with opposite spins are spatially separated on alternating 
A and B sites. In other words, it represents how much 
those electrons respect the staggered order. Upon the 
same gauge transformation with matrix (C4), 



R a 



= & H a. (C7) 



The SU(2) gauge field associated with Eq. (Cl) is bro- 
ken into an Abelian magnetic field due to effective Higgs 
mechanism, 

Fij = dii&.Aj) - dj(<5>-Ai) + 2$ • (ft* x 0,-*), (C8) 



B - -e -bF b - — 



(C9) 
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where <I> = &h/\&h\, and Eq. (C9) is the magnetic field 
of a Dirac monopole. 

An unsolved issue is the Bogomol'nyi relation 18,23 . For 
the non-Abelian gauge field f^j = diAj —djAi — i[Ai, Aj], 
and the covariant derivative T>i = di — i[Ai, ], it is 
straightforward to derive 



= ^ijk[Dk<i>H - (| - | + J?2)3/2 <Tk]. (cio) 

If not were the last term, Eq. (CIO) reproduces the Bo- 
gomol'nyi relation. One can show that only for a profile 
function = 2R/ sinh(i?) (the case of true 't Hooft- 
Polyakov monopole) that the last term vanishes. While 
our £(R) asymptotically resembles 2R/ sinh(,R), it gives 
a different profile for finite R. As a result, the last term 
in Eq. (CIO) vanishes only when R — >• oo. 

Appendix D 

We eliminate dt on both sides of Eq. (8a), 

ds = (1 -£ 2 )(s -n)dn. (Dl) 



Firstly, dot product both sides of Eq. (Dl) with n, re- 
garding the constraint n 2 = 1, 

n • ds = (1 - £ 2 )(s • n)^d(n 2 ) = 0, (D2) 

which gives us the relation d(n • s) = s • dn. Secondly, 
dot product with s on both sides of Eq. (Dl) we get, 

d(s 2 ) = 2(1 - £ 2 )(s • n)(s • dn) = (1 - £ 2 )d(s • n) 2 

(D3) 

From Pythagorean theorem we know that s 2 = (s • n) 2 + 
(s x n) 2 , then take derivative on both sides, regarding 
Eq. (D3), we arrive at 

£ 2 d(s-n) 2 = -d(s x n) 2 . (D4) 

Assume the initial condition (s ■ n)\o = 1, the above 
equation can be integrated into, 

(s-nf+ { -^pl = sl + ^±^ = l, (D5) 

which justifies Eq. (10). 
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